BAB4R-CONF-03 /002 
SLAC-PUB-9667 
March 2003 



Rates, Polarizations, and Asymmetries in Charmless 
Vector- Vector B Decays 



The BABAR Collaboration 
March 13, 2003 

Abstract 

With a sample of approximately 89 million BB pairs collected with the BABAR detector, we measure 
branching fractions, determine the degree of longitudinal polarization, and search for direct CP 
violation in the decays B° — ► cj)K* and B + — > <pK* + . We perform a search for other charmless 
vector- vector B decays involving p and if* (892) resonances and observe the decays B + — > p°K* + 
and B + — > p°p + . The branching fractions are measured to be B(4>K*°) = (11.1 jli'2 ^ x 
10" 6 , B(<f>K*+) = (12.1^-1 ± L5 ) x 10 -6 j B{p°K* + ) = (7.7±|;5 ± 1.4) x lO" 6 , and B(p°p + ) = 
(9.9^2'5 =t 2.5) x 10 -6 . The longitudinal polarization fractions are measured to be Tl/T((()K* ) = 
0.65 ± 0.07 ± 0.04 and T L /T((t)K* + ) = 0.46 ± 0.12 ± 0.05. We measure the charge asymmetries: 
A C p{4>K*°) = +0.04 ± 0.12 ± 0.02 and A C p{(f>K* + ) = +0.16 ± 0.17 ± 0.04. 
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1 INTRODUCTION 



The first evidence for charmless vector-vector B decays was provided by the CLEO experiment with 
the measurement of one B — > <j)K* channel [Q. Measurement of both B —* cj)K* charge modes was 
performed by BABAR 2 , which also reported the search for direct CP violation in these modes 
The CLEO experiment also set upper limits on the decay rates in several other vector-vector final 
states of B decays [I]. The Belle experiment announced large signal yield in the B + — ► p°p + 
channel 

Recently, there has been interest in charmless B decays because of the clean environment for 
the search for new physics. For example, new particles contributing to penguin diagrams, such 
as charged Higgs bosons or SUSY particles, would provide additional amplitudes with different 
phases. Charmless B decays are also sensitive to the weak phases a = arg [ — V td V t l / V ud V* b ] and 
7 = arg [ — V ud V* b / V cd V* b ] arising from the elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix [0]. The vector- vector charmless B decays provide additional information about 
the decay dynamics and strong phases, which could be obtained from the analysis of angular 
distributions jjj. 

The decays B — > (pK* are expected to proceed through pure penguin diagrams (b — > s loops) as 
illustrated in Fig.^ Similarly, the decays B — ► pK* are expected to be dominated by b — > s penguin 
transitions with additional contributions from Cabibbo-suppressed tree-level b — > u transitions, 
while the decays B —* pp proceed primarily through Cabibbo-favored tree-level b — > u transitions 
and CKM-suppressed b — > d penguins. 



B b 




u,d u,d u,d u,d 



Figure 1: Two of the dominant diagrams describing the decays B — ► pp, pK* , and <j>K* . 



The measurement of direct CP violation in pure penguin modes, such as B — > 4>K* , is sensitive 
to non-standard-model physics jHj. In the standard model, direct CP violation could arise due to 
the difference between the b — > u tree and b — > s (b — * d) penguin amplitude weak phases (Hj, which 
is 7 (a) in the case of the decays B — > pK* (B — ► pp). Direct CP violation would be observed as 
an asymmetry of B decay rates: 

= T(B^f)-T(B^f) 

CP r(B-/) + r(s-/) ' {) 

However, large uncertainties in the strong phases, which can be calculated by certain models, 
weakens the quantitative relationship to the weak phases. 
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The time-dependent asymmetries in B decays to CP eigenstates would provide important tests 
of the standard model ^Uj. Comparison of the value of sin 2/3 obtained from (pK*°(^ K^ir ) with 
that from charmonium modes can probe for new physics. Time-dependent measurements in B — > pp 
modes combined with isospin relations among the decay amplitudes for these modes would provide 
a measurement of a. Angular analysis is important for time-dependent asymmetries because of the 
mixture of CP-odd and CP-even components, and for the isospin analysis of B — > pp modes. 

2 THE BABAR DETECTOR AND DATASET 

In this analysis we use the data collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ collider [2] located at the Stanford Linear Accelerator Center. The results presented 
in this paper are based on data taken in the 1999-2002 run comprising an integrated luminosity 
of 81.9 fb _1 , corresponding to 88.9 million BB pairs, at the T{AS) resonance (on-resonance) and 
9.6 fb _1 approximately 40 MeV below this energy (off-resonance). The T(AS) resonance occurs 
at the e + e~ center-of-mass (cm.) energy, y/s, of 10.58 GeV. The asymmetric beam configuration 
in the laboratory frame provides a boost to the T(4S) increasing the momentum range of the B 
meson decay products up to about 4 GeV/c. 

Charged-particle momenta are measured in a tracking system that is a combination of a silicon 
vertex tracker (SVT) consisting of five double-sided detectors and a 40-layer central drift chamber 
(DCH), both operating in a 1.5-T solenoidal magnetic field. BABAR achieves an impact parameter 
resolution of about 40 pm for the high momentum charged particles in the B decay final states, 
allowing the precise determination of decay vertices. The tracking system covers 92% of the solid 
angle in the cm. frame. The average track-finding efficiency is approximately 98% for momenta 
above 0.2 GeV/c when the angle between track and the beam axis is greater than 500 mrad. 

Charged particle identification is provided by the average energy loss (dE/dx) in the tracking 
devices (SVT and DCH) and by an internally reflecting ring imaging Cherenkov detector (DIRC) 
covering the central region. A K-tt separation of better than four standard deviations (a) is 
achieved for momenta below 3 GeV/c, decreasing to 2.5<r at the highest momenta in the B decay 
final states. Photons are detected by a CsI(Tl) electromagnetic calorimeter (EMC), which provides 
excellent angular and energy resolution with high efficiency for energies above 30 MeV. For a 1 GeV 
photon the energy and angular resolutions are 3% and 4 mrad, respectively. Electrons are identified 
by the tracking system and the EMC. 

3 EVENT SELECTION 

Hadronic events are selected based on track multiplicity and event topology. We fully reconstruct 
B meson candidates from their charged and neutral decay products including the intermediate 
states 4> -> K+K~, K*° -> K+tt', K*° -> K°7r°, K*+ -> K + ir°, K*+ -> K°vr+, p° -> Tr+vr", 
p + — > 7T + 7r°, with 7T° — ► 77 and K° — ► Kg — ► tt + tt~ ^H]- Candidate charged tracks are required to 
originate from the interaction point, and to have at least 12 DCH hits and a minimum transverse 
momentum of 0.1 GeV/c. Looser criteria are applied to tracks forming K s candidates to allow for 
displaced decay vertices. The K s candidates are required to satisfy \m n + n - — rn K o\ < 12 MeV/c 2 
with the cosine of the angle between their reconstructed flight and momentum directions greater 
than 0.995 and the measured proper decay time greater than five times its uncertainty. Kaon tracks 
are distinguished from pion and proton tracks via a likelihood ratio that includes, for momenta 
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below 0.7 GeV/c, dE/dx information from the SVT and DCH, and, for higher momenta, the 
Cherenkov angle and number of photons as measured by the DIRC. 

We reconstruct ir° mesons from pairs of photons, each with a minimum energy 30 MeV. The 
typical width of the reconstructed 7r° mass is 7 MeV/c 2 . A ±15 MeV/c 2 interval is applied to select 
7T° candidates. We select (ft, K*, and p candidates with the following requirements on the invariant 
masses of their final states: 0.99 < m K + K - < 1.05, 0.75 < rriKw < 1-05, and 0.52 < < 1.00 (all 
in GeV/c 2 ). The helicity angle # x of a (j), K*, or p (x = 1, 2 for the two resonances in the B decay) is 
defined as the angle between the direction of one of the two daughters {K + , K, or tt + respectively) 
and the parent B direction in the resonance rest frame. To suppress combinatorial background we 
restrict the K* — > Kir and p + — ► 7r + 7r° helicity angles (cos# x > —0.5). This effectively requires 
the 7T° momentum to be larger than 0.35 GeV/c. 

We identify B meson candidates kinematically using two nearly independent variables |llj . 
the beam energy-substituted mass toes = [(s/2 + pi ■ ps) 2 / Ef — p^] 1 ^ 2 and the energy difference 
AE = (EiEs — Pi ■ Pb — s /2)/\/s, where (Ei, pi) is the initial state four-momentum obtained from 
the beam momenta, and (Eb,Pb) is the four-momentum of the reconstructed B candidate. For 
signal events toes peaks at the B mass and AE at zero. Our initial selection requires toes > 5.2 
GeV/c 2 and \AE\ < 0.2 GeV. 

Charmless hadronic modes suffer from a large background due to random combinations of tracks 
produced in quark-antiquark continuum events (e + e~ — > qq, q = u,d, s, c). Background events from 
the continuum are distinguished by a jet-like structure as opposed to the more spherical topology 
of BB pairs produced in T(4S) events. To reject continuum background we require | cos#t| < 0.8, 
where 0t is the angle between the thrust axis of the B candidate and that of the rest of the 
tracks and neutral clusters in the event, calculated in the cm. frame. The distribution of the 
cos 6t variable is sharply peaked near ±1 for combinations drawn from jet-like qq pairs, and nearly 
uniform for the isotropic B meson decays. We also construct a Fisher discriminant that combines 
eleven variables 14 : the polar angles of the B momentum vector and the £>-candidate thrust axis 
with respect to the beam axis in the T(4S) frame, and the scalar sum of the cm. momenta of 
charged particles and photons (excluding particles from the B candidate) entering nine coaxial 
angular intervals of 10° around the 5-candidate thrust axis. 

Monte Carlo (MC) simulation JH] demonstrates that contamination from other B decays is 
negligible for the modes with narrow (ft resonance and is relatively small for other charmless B decay 
modes. We achieve further suppression of B decay background by removing all signal candidates 
that have decay products consistent with D — > Kir, Ktttt decays. The remaining small background 
coming from B decays is accounted for in the fit. It is found that background subtraction is 
necessary only in the B — » pK* analysis, where there is a small probability of charmed B decays 
being reconstructed as signal. In this analysis we assume negligible contribution of other partial 
waves in our final states selected within vector resonance mass windows. 



4 ANALYSIS METHOD 

We use an unbinned extended maximum likelihood (ML) fit to extract signal yields, charge asymme- 
tries, and angular polarizations simultaneously. We define the likelihood for each event candidate: 

3 2 

£i = ^2^2n jk V jk (xi;a), (2) 
j=i k=i 
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where Vj^ixi, a) is the probability density function (p.d.f.) for measured variables x*j of a candidate 
i in category j and flavor state k, and rijk are the yields to be extracted from the fit. There are 
three categories: signal (j = 1), continuum qq (j = 2), and BB combinatorial background (j = 3). 
The p.d.f. 's are non-zero only for the correct final state flavor (k = 1 for B —* / and k = 2 for 
-B — » /). The i? flavor is determined by its charge, except for the 4>K*° final state where the 
flavor is determined from the charge of the kaon from the K*° — ► ET + 7r - decay, and the flavor is 
not defined for the K*° — > i^°7r° decay. The fixed parameters a define the expected distributions 
of measured variables in each category and flavor state. We rewrite the event yields njk in each 
category in terms of the asymmetry Aj and the total event yield ny. nj\ = rij x (1 + Aj)/2 and 
rij2 = n,j x (1 — Aj)/2. This definition is consistent with Eq. ^ 

The fit input variables a?j are AE, m-ES) Fisher discriminant, invariant masses of the K* and (j) 
(or p) resonances, and the K* and <f> (or p) helicity angles # x (x = 1, 2). The correlations among the 
fit input variables in the data and signal MC are found to be small (typically less than 5%). The 
p.d.f. Vjk(xi] a) for a given candidate i is the product of the p.d.f.'s for each of the variables, except 
for the helicity angles. We take into account the angular correlations in the signal and the detector 
acceptance effects in the helicity angle p.d.f. parameterization. Due to the limited statistics in 
our analysis we adopt a simplified angular analysis technique where we integrate over the angle 
between the decay planes of the two vector-particle decays, leaving a p.d.f. that depends only on 
the two helicity angles. This distribution is sufficient to determine the longitudinal polarization 
fraction fi = Tl/T. The differential decay width is defined as: 

h fj 7, =^(j(l-/L)sin 2 g 1 sin 2 g 2 + / L cos 2 g 1 cos 2 g 2 j . (3) 

T a cos 6 \d cos 02 4 L4 J 

We allow for multiple candidates in a given event by assigning candidates a weight of 1/iVj, 
where iVj is the number of candidates in the same event. The average number of candidates per 
event is close to one (varying from 1.05 to 1.14 depending on the mode). The extended likelihood 
for a sample of iV can d candidates is 

eH 1 ^)- < 4 » 

The event yields rij, asymmetries Aj, and polarization fi are obtained by minimizing the 
quantity x 2 = ~ 21n£ using the minimization package MINUIT |16j . The dependence of x 2 on a fit 
parameter rij, Aj, or Jl is obtained with the other fit parameters floating. We quote statistical 
errors corresponding to unit change in \ 2 . When more than one K* decay channel is measured 
for the same B decay, the channels are combined by adding their x 2 distributions for rij, Aj, or 
//,. The statistical significance of a signal is defined as the square root of the change in x 2 when 
constraining the number of signal events to zero in the likelihood fit; it describes the probability 
for the background to fluctuate to the observed event yield. 

The fixed parameters a describing the p.d.f.'s for signal and background distributions are ex- 
tracted from MC simulation, on-resonance AE-ities sidebands, and off-resonance data. The MC 
resolutions are adjusted by comparisons of data and simulation in abundant calibration channels 
with similar kinematics and topology, such as B — > Dn, Dp with D — ► Ktttt, Kit. The simulation 
reproduces the event-shape variable distributions found in data. 

To describe the signal distributions, we employ Gaussian functions for the parameterization 
of the p.d.f.'s for AE and ttjes an d a relativistic P-wave Breit-Wigner distribution convoluted 
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with a Gaussian resolution function for the resonance masses. For the background we use low- 
degree polynomials or, in the case of toes) an empirical phase-space function ^7j. The background 
parameterizations for resonance masses also include a resonant component to account for resonance 
production in the continuum. The background helicity angle distribution shape is again separated 
into contributions from combinatorics and from real mesons, both fit by low-degree polynomials 
multiplied by an empirical function 1/(1 + exp((# x — 0q)/o)) to account for the detector acceptance 
effects, where 9q and a are fixed parameters. For both the signal and background, the p.d.f. for the 
Fisher discriminant is represented by a Gaussian with different widths above and below the mean. 

5 PHYSICS RESULTS AND SYSTEMATIC UNCERTAINTIES 

The results of our maximum likelihood fits are summarized in Tabled For the branching fractions 
we assume equal production rates of B°B° and B + B~~ . We find significant signals in both cf)K* 
decay modes. We measure the charge asymmetries and longitudinal polarizations in all <f)K * final 
states. We also observe a significant yield of events in p °K* + (4.7a) and p°p + (4.4a) final states. 
The projections of the fit results are shown in Fig. |2] and |31 where we plot only a subsample of 
events, enhancing the signal with a requirement on the ratio of the signal probability to background 
probability ("P s ig and "Pbkg from Eq. [2J) . 

We study the performance of the ML fit with the MC samples where the signal events are taken 
from the complete MC simulation and the background is distributed according to the sideband 
parameterizations. The input values for the number of signal events, charge asymmetry, and 

Table 1: Summary of results for the measured B decay modes; e denotes the reconstruction effi- 
ciency and etot the total efficiency including daughter branching fractions, both in percent; n s i g is 
the fitted number of signal events, B is the branching fraction, Acp is the signal charge asymmetry, 
and Jl is the longitudinal polarization. The decay channels of K* are shown when more than one 
final state is measured for the same B decay mode. All results include systematic errors, which are 
quoted following the statistical errors. 



Mode 


e 


Etot 




B (xHT 6 ) 


Acp 


h 






5.3 




12.1 tf l ± 1.5 


+0.16 ±0.17 ±0.04 


0.46 + 0.12 + 0.05 




26.0 
14.3 


2.9 
2.3 


33.3lg 4 ± 1.2 
22.3±e'.5 ± 3.2 


12.7 tH ± 1-2 
10.7 m ± L9 


-0.02 ±0.20 ±0.03 
±0.63 to.ll ±0.05 


0.50 ±015 ± 0.04 
0.40 ± 0.07 


<f>K*° 




10.2 




11.1 ±{\l ± 1.1 


+0.04 ±0.12 ±0.02 


0.65 ±0.07 ±0.04 




29.7 
11.4 


9.7 
0.6 


101+ " ±3 
2.0^3 ± 0.6 


11.7 ±1.4 ±1.1 
3.5 t^il.l 


+0.04 ±0.12 ±0.02 


0.64 ±0.07 ±0.03 
1.00 ±°-gg ± 0.25 


p°K*+ 




8.4 




7.711:1 ± 1-4 








21.0 
10.9 


4.8 
3.6 


44.4+1^ ±3.4 
9.1±JV ±5.2 


10.4+2;? ± 1.7 
2.9tg.o ± I- 8 






P°P + 


11.3 


11.1 


97.5la4.3 ± 12 -! 


9.912.5 ± 2 - 5 
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5.200 5.225 5.250 - 5.275 5.300 5.200 5.225 5.250 - 5.275 5.300 

m FS (GeV/c ) m FS (GeV/c ) 



Figure 2: Projections onto the variable m ES for B° -> (j)K*° (a), B + -> (b), -> p°i-C*+ 

(c), and B + — > ("dj candidates after a requirement on the signal-to-background probability 

ratio 'Psig/'Pbkg with the p.d.f. for tties excluded. The histograms show the data, which are the sum 
of the two K* decay channels when appropriate, while the shaded area is K* — > Kir channel alone. 
The solid (dashed) line shows the signal-plus-background (background only) p.d.f. projection. 
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Figure 3: Vector resonance invariant mass distributions for B + — ► p°K* + candidates: m n + n - (a) 
and rriK-K (b); and for B + — > p°p + candidates: m^+ n - (c) and m n + n o (d), after a requirement 
on the signal-to-background probability ratio 'Psig/'^bkg with the p.d.f. for the invariant mass 
excluded. The histograms show the data and the shaded area corresponds to the p°K* + final 
state with K* + — ► K + tt° channel alone. The solid (dashed) line shows the signal-plus-background 
(background only) p.d.f. projection. 
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decay polarization are well reproduced, and any small discrepancy is accounted for in the systematic 
errors. 

Systematic uncertainties in the ML fit originate from assumptions about the signal and back- 
ground distributions. Uncertainties in the p.d.f. parameters arise from the limited statistics in 
the background sideband data and signal control samples. We vary the p.d.f. parameters within 
their respective uncertainties, and derive the associated systematic errors. The signals remain sta- 
tistically significant under these variations. Additional systematic errors in the number of signal 
events originate from uncertainty in the charm background subtraction in the pK* channels where 
we take the uncertainty to be 100% of the subtracted values (2 and 5 events in K* + — » K°tt + and 
K* + — > K + ir° channels respectively). 

The dominant systematic errors in the efficiency are track finding (0.8% per track), particle 
identification (2% per track), and K$ and tt° reconstruction (5% each). Other minor systematic 
effects from event selection criteria, daughter branching fractions JHj> MC statistics, and B meson 
counting sum to less than 4%. The efficiency in the ML fit to signal samples can be less than 
100% because of fake combinations passing the selection criteria, and we account for this with a 
systematic uncertainty (3-12%). This effect is larger in the final states with p resonances because 
of the broader distributions. Efficiency uncertainties affect the value of the branching fraction, but 
not its significance. 

In the search for pK* and pp final states we fit only for event yields and exclude angular and 
S-flavor information from the analysis. The reconstruction efficiency depends on the decay polariza- 
tion. We average reconstruction efficiencies for the 100% transverse and 100% longitudinal angular 
polarizations for each decay channel and assign the systematic errors as the root-mean-square of 
the uniform efficiency distributions between the two extreme cases. The resulting systematic errors 
are 9% in p°K* + (-^ K°tt + ), 19% in p°K* + (^ K + tt°), and 18% in p°p + decay channels. For the 
(pK* final states we calculate the efficiencies assuming an average polarization of (0.60 ±0.06 ±0.05) 
and assign a systematic error corresponding to the total polarization measurement error. 

In the polarization measurements, we again include uncertainties from p.d.f. variations where 
we account for possible bias in the detector acceptance and background parameterizations. The 
biases from the finite resolution in helicity angle measurement and dilution due to the presence 
of the fake combinations are studied with MC simulation and are accounted for with conservative 
systematic error of 0.03. 

We find the charge asymmetry of the track reconstruction efficiency to be consistent with 
zero within an uncertainty of less than 0.01 for a wide range of momenta [3;. Taking into account 
particle identification requirements similar to the ones applied to the K* daughters, the asymmetry 
is consistent with zero with an uncertainty of 0.02. We also find a negligible effect on the measured 
asymmetry from any possible bias in track-momentum measurements studied in e + e~ — > p + p~ and 
cosmic ray events. The asymmetry measurement in the B° — ► <f>K*° decay mode is corrected by the 
inverse dilution factor 1/(1 — 2w), where w, the fraction of doubly misidentified Kir combinations 
originating from K*°, is less than 0.01. 

6 SUMMARY 

We have measured branching fractions, longitudinal polarizations, and charge asymmetries with 
the decays B° — ► (j)K*° and B + — > <j)K* + . Because the final states contain three strange quarks or 
antiquarks, in the standard model they are necessarily due to penguin diagrams. This makes them 
particularly susceptible to non-standard-model contributions. We observe the decays B + — * p°K* + 
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and B + —* p°p + and report the corresponding branching fractions. The B + — > p°K* + decay process 
is of interest for direct CP measurements due to possible large penguin-tree interference, while the 
B — > pp decays have potential for the measurement of the weak phase a. 

The B — * (j)K* branching fractions are in agreement with our earlier less precise measure- 
ments 0. Our charge asymmetry results rule out a significant part of the physical Aop region, 
allowing for constraints on new physics models jH], but are not yet of sufficient precision to allow 
precise comparison with standard model predictions 19.. We have performed the angular analysis 
in the penguin-dominated rare B decays and measure a relatively large longitudinal polarization 
in the decay amplitude, as predicted 7,. Our measurement of B + — ► p°p + branching fraction is 
significantly lower than the central value measured by Belle (Hj even taking into account polariza- 
tion uncertainty. Our results are preliminary and use increased statistics compared to the earlier 
BABAR measurements of the B — * <f)K* branching fractions and charge asymmetries [2 |Hj . 
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